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Abstract 

Thin films are usually obtained by depositing atoms with a continuous flux. 
We show that using a chopped flux changes the growth and the morphology of 
the film. A simple scaling analysis predicts how the island densities change as 
a function of the frequency of the chopped flux in simple cases where aggrega- 
tion is irreversible. These predictions are confirmed by computer simulations. 
We show that the model can be used to obtain information on the diffusion 
or the evaporation of the adatoms. The model is also useful to understand 
the growth of thin films prepared by pulsed sources. 
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One of the main interests of usual deposition techniques such as Molecular Beam Epitaxy 
|jlj] is that the structure of the deposited films is to a large extent determined by kinetic 
factors, as opposed to thermodynamic equilibrium. This allows to "play games" 0] with the 
different growth parameters (incident flux of particles, diffusion coefficient of an adatom . . . ) 
in order to obtain different film morphologies. A simple example is given by the quantity 
of islands grown on a substrate at low enough temperatures : it is known that the number 
of islands at saturation is given by (F/D) 1 ^ 3 |3|-[5| where F is the incident flux and D the 
diffusion coefficient. Then, by increasing the flux or decreasing the diffusion constant (by 
lowering the substrate temperature), one can adjust the saturation number of islands grown 
on the substrate. In this sense, each kinetic factor is a "handle" on the system, allowing to 
control the morphology of the films. We introduce in this Letter a new kinetic handle, which 
should enable a larger control over film growth : the chopping of the incident flux. We note 
that this flux modulation is intrinsic to other deposition techniques such as cluster laser 
vaporization (the laser is pulsed ||). It is therefore important to understand how growth 
proceeds in the presence of a modulated flux if one is to be able to interpret experiments 
performed in these conditions. For example, one may wonder whether the usual growth 
theories |j|-[5| can be used by replacing the continuous flux by the average value of the 
chopped flux over a cycle. In the following, we will show that this is not the case, and that 
the growth of the film is profoundly changed by the modulation of the incident flux for the 
case of growth with irreversible aggregation (critical island size 1, see [^H^J^HH- Conversely, 
we show what kind of information can be derived from experiments carried under these 
conditions. 

The basic idea of our method is that if instead of using a continuous flux we use a 
chopped flux to grow a film, the number of islands formed on a substrate will depend on the 
chopping frequency / and on d, the fraction of the period the flux is "on" (see Fig. [[]). 
This dependence is due to the fact that the free particle concentration on the surface does 
not reach its steady state concentration instantaneously, but only after a time which we will 
call r m . Then, if the timescale of the chopping (1/f) is much smaller than r m , the system 
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only sees the average flux. In the contrary case, everything happens as if the instantaneous 
flux was used instead. Then, there will be a transition from one behaviour to the other at 
a chopping frequency close to l/r m . We discuss this in detail in next section in the case of 
growth without desorption: the corresponding analysis in presence of evaporation is briefly 
presented afterwards. 

Presentation of the method : Fig |T] shows schematically the time evolution of the free 
adatom (monomer) concentration in presence of a modulated flux. We assume here, for 
simplicity, that adatom-adatom aggregation is irreversible and that only monomers can 
move on the surface f|-|5],[7|-§| . Then, the important timescale for growth is given by the 
time needed for monomers to reach their steady state concentration or to disappear. Since, 
in the absence of evaporation, the monomers disappear mainly by diffusing randomly until 
they aggregate with an island ||, their mean lifetime on the surface is given by r m ~ l 2 /D 
where 21 is the mean distance between islands and D the diffusion constant of the monomers. 
We obtain r m ~ 1/ (DN) where N is the island concentration. We can predict three regimes 
of behaviour, depending on the relative magnitude of r m , d/f and 1/f. 

For low chopping frequencies (r m <C d/f), the monomer concentration reaches its steady 
state concentration p ss = Fr m almost instantaneously in the timescale of a period. After 
the flux is turned off, the monomer concentration goes back to also almost instantaneously 
ijm "C 1//)- Then, between two successive "flux on" periods, nothing happens since only 
the monomers can move, and there is no monomer left. Therefore, growth proceeds as if we 
had a continuous flux Fi and the island concentration at saturation for low frequencies Nif 
satisfies the well-known result ~ (Fi/D) x , with \ ~ 0.36 for fractal islands 0. 

When the frequency increases (d/f <C r m ), the monomer concentration cannot attain its 
steady state concentration during the "flux on" part of the cycle : pit = dj f) ~ Fid/ f <C p ss . 
Since we are interested in d <C 1, p vanishes before next cycle, and the average of 
the monomer concentration over a cycle is smaller than in the preceding regime, lead- 
ing to a smaller island density. More precisely, the increase of island density during 
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a cycle is given by : AiV = J cyde Dp 2 (t)dt. Using p(t) = (Fid/ f)e */ Tm , we obtain 
AN ~ D J^ f (Fid/f)e- 2t ^ ~ 1/2 (F^/ ff * (Dr m ). Since r m ~ we obtain 

A(7V 2 ) ~ (Fid/ff. We have checked that, as in the case of a continuous flux PJlQl, the 



saturation island density is reached for a surface coverage of about 15%. Then, the total 
number of cycles needed to reach saturation is n c ~ //(6i^d). Since A(iV 2 ) is a constant 
over the cycles, we simply multiply A(N 2 ) by the number of cycles to obtain iVj the satu- 
ration island density in this intermediate regime : Ni ~ (Fid/f) 1/2 . We note that here the 
saturation island density does not depend on the diffusion constant. This regime is limited 
by the fact that if r m ~ 1//, the monomer concentration no longer fully decays when the 
next "on" cycle begins. The monomer concentration slowly increases, reaching its steady 
state value after several cycles. 

We then cross to the regime of high frequencies (r m ^> l/f), where many deposition 
cycles are carried out during the monomer equilibration, and the system only sees the 
average flux F av . Then the island concentration at saturation for high frequencies N^f 
satisfies N hf ~ (F av /D) x = d x N lf < N lf . 

Measuring the evaporation time : We briefly treat the case of film growth in presence of 
significant evaporation of monomers from the surface (for a more detailed treatment in the 
presence of a continuous flux, see fL0|). In this case, the lifetime of the adatoms is no longer 
determined by the aggregation with islands but rather by random evaporation after a mean 
time r e PJ5|flO[1. We study here the dependence of the island density N e at a fixed amount of 
deposited matter as a function of the modulation frequency. This quantity is not necessarily 
the saturation island density but it is easier to measure experimentally when evaporation is 
important. As in the non evaporation case, three different regimes will be scanned as the 
frequency increases. 

For low chopping frequencies (r e d/f), the monomer concentration reaches its steady 
state concentration p ss = Fr e and vanishes almost instantaneously in the timescale of a 
period. Therefore, during a cycle, the island density increases according to dN e /dt ~ Dp 2 ~ 
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DFfr^, leading to AN e (per cycle) = DF^r^d/f. The island density for fixed M = F av t in 
this regime is therefore Nff ~ DFiT e M since the number of cycles is n c = fMf(dFj). 

When the frequency increases (d/f <C r e ), the monomer concentration cannot attain its 
steady state concentration during the "flux on" part of the cycle : pit = dj f) ~ Fid/ f <C p ss - 
Again, the increase of island density during a cycle is given by : AN e = f cycle Dp 2 (t)dt ~ 
1/2 (Fid/f) 2 Dr e . Counting the number of cycles as before leads to an island density in 
this intermediate regime iVf ~ DFidr e M f~ l . As in the case of complete condensation, this 
regime is limited by the fact that if r e ~ 1/f, the monomer concentration no longer fully 
decays when the next "on" cycle begins. Then, the monomer concentration slowly increases, 
reaching its steady state value after several cycles. 

In the regime of high frequencies (r e ^> 1/f), many deposition cycles are carried out 
during the monomer equilibration, and the system only sees the average flux F av . Using the 
result obtained in the low frequency case (see also Ref. |10j), the island concentration at 
saturation for high frequencies can be written as ~ DF av r e M = dNfj <C Nf*. 

Computer simulations : Monte-Carlo simulations of the growth of thin films by a con- 
tinuous flux have been independently implemented by several groups with small variations 



TT| |13|| . To simulate the growth in the presence of a chopped flux, we introduce the four 
following physical processes : 

(1) Deposition. Particles are deposited at randomly-chosen positions of the surface at a 
flux Fi during the "on" fraction of a cycle (d <C 1). During the rest of the cycle, no particle 
reaches the surface (see Fig [l]). The average flux reaching the surface is thus F av = dFi. 

(2) Diffusion. Each diffusion time r, all monomers are chosen at random and moved in 
a random direction by one lattice site. 

(3) Evaporation. Isolated adatoms can evaporate off the surface at a constant rate. We 
denote by r e the mean lifetime of a free adatom on the surface. 

(4) Aggregation. If two particles come to occupy neighboring sites, they stick irreversibly 
and form an island. Islands are assumed to be immobile. 
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We first analyse the results when desorption is negligible (r e = oo). Fig |2| shows the 
evolution of the island density for a coverage of 15% (which we have checked is also the 
maximum island density, see also 0j) as a function of the rescaled frequency f/d. We see 
the three regimes predicted by our simple analysis. In particular, the relation iVj ~ f~ 1 ^ 2 is 
well verified in the intermediate regime. The inset shows that the ratio of the low and high 
frequencies island densities scales as d x , as predicted. 

When desorption is important, the agreement is also very good. Fig ||| shows that the 
relation Nf ~ f~ x is well verified in the intermediate regime. As in the case of complete 
condensation, the inset in Fig |3| shows that the ratio of island densities in the low and high 
frequencies regimes scales as d. 

Summary, Perspectives : As anticipated, we have shown that it is not possible to interpret 
growth experiments under a chopped flux in the framework of usual growth models 10-0] 
which assume a continuous flux. For example, the relation of the saturation island density 
to the experimental parameters (flux, diffusion coefficient) depends on the regime, i.e. on 
the frequency of the chopped flux. 

What kind of information can be derived from flux-modulated experiments for usual 
deposition techniques such as Molecular Beam Epitaxy? We note that modulation of the 
flux can be achieved rather easily by intercalating a rotating disk with a slit close to the 
sample. The order of magnitude of the frequencies required to observe the three regimes is 
reasonable. For example, in Stroscio's Fe/Fe(001) experiments [14|, at room temperature, 



N ~ 7 10 12 cm" 2 and D ~ 7 10- 12 cm 2 s _1 then f high ~ ADN ~ 2mHz and f low = df h ~ 
1Hz for d = 0.01. In the case of complete condensation, measuring the saturation island 
density N sat as a function of frequency can help checking the growth conditions. First, 
N sat should show a sharp decrease for frequencies close to l/r m . Second, the ratio of the 
densities in the two extreme regimes should scale as d x if aggregation is indeed irreversible. 
We notice that this test - which is similar to the N sat versus flux test [f^|?|,|8|,|T5|| - can be 
carried out experimentally without changing the flux, or even without knowing it. The 
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effects of reversible aggregation should be addressed in future studies. The intermediate 
regime is also interesting since N sat is independent of the diffusion constant (and therefore 
of temperature). In growth experiments carried under a continuous flux, this constant 
density may only be observed for growth in presence of trapping impurities or reconstructed 
surfaces. In the case of growth in presence of important desorption, we speculate that this 
method allows to measure the evaporation time in-situ. By using reflection-high-energy- 
electron diffraction (RHEED) one could in principle measure a characteristic time needed for 
growing a monolayer, say a quarter oscillation. Although RHEED oscillations are difficult to 
interpret in detail, we anticipate that this characteristic time should depend strongly on the 
modulation frequency in the neighborhood of / ~ l/r e . Therefore, studying the evolution 
of this characteristic time as a function of the chopping frequency may lead to an estimate 
of T e . 
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FIGURES 

FIG. 1. Schematic illustration of the chopped flux (bold) and the corresponding evolution of the 
monomer density (dashed line) . F, and F av refer to the instantaneous and average flux respectively. 

FIG. 2. Monte Carlo simulation of the island densities as a function of the rescaled chopping 
frequency f/d in the case of complete condensation. The island densities correspond to a surface 
coverage of 15%, an instantaneous flux Fj = 10~ 7 and a diffusion constant D = 1/4. Each curve 
corresponds to a different value of d : d = .5 (circles), d = .1 (squares), d = .01 (diamonds), 
d = .001 (triangles). The solid line has a slope of -1/2. The inset shows the dependence of the high 
frequency island density on d. The solid line is a fit to the data and has a slope 0.36, in excellent 
agreement with the predicted slope (% = 0.36). 

FIG. 3. Computer simulation of the island densities as a function of the rescaled chopping 
frequency f/d when growth is limited by evaporation. The island densities correspond to a depo- 
sition time F av t = 0.77, an instantaneous flux Fj = 10 -7 , a diffusion constant D = 1/4 and an 
evaporation time r e = 100. Each curve corresponds to a different value of d d = .5 (left triangles), 
d = .1 (squares), d = .01 (circles), d = .003 (triangles up), d = .001 (diamonds). The solid line has 
a slope -1. The inset shows the dependence of the high frequency island density on d. The solid 
line is a fit to the data and has a slope 0.97, in excellent agreement with the predicted slope of 1. 
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